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ABSTRACT 

o" 

^^ I We report the discovery of a new double-image gravitational lens resulting from our 

search for lenses in the southern sky. Radio source PMN J2004-1349 is composed of 
two compact components separated by 1'.'13 in VLA, MERLIN and VLBA images. The 

Q ■ components have a flux ratio of 1:1 at radio frequencies ranging from 5 GHz to 22 GHz. 

The /-band optical counterpart is also an equal double, with roughly the same separation 
and position angle as the radio double. Upon subtraction of the components from the 

pg , /-band image, we identify a dim pattern of residuals as the lens galaxy. While the 

present observations are sufficient to establish that PMN J2004-1349 is a gravitational 

I/-) | lens, additional information will be necessary (such as the redshifts of the galaxy and 

quasar, and precise astrometry and photometry of the lens galaxy) before constructing 

q ■ detailed mass models. 
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1. Introduction 

The new double-image gravitational lens presented in this paper was found during a systematic 
search for lenses in the southern sky. Among the reasons to search for new lenses are: to find 
objects well-suited for determining the Hubble constant via time-delay measurements (Refsdal 1964; 
Myers 1999; Schechter 2000); to constrain the cosmological constant by measuring the incidence of 
multiple- image lensing (Turner 1990; Fukugita et al. 1992; Kochanek 1996; Falco et al. 1998); and 
to study the mass distributions (e.g., Keeton, Kochanek & Seljak 1997; Rusin & Ma 2000) and 
evolution (Keeton, Kochanek & Falco 1998) of medium-redshift galaxies. 

Our search methodology will be described in a future paper, but is briefly as follows. Our 
sample consists of southern radio sources (0° > 5 > —40°) with flat radio spectra (a > —0.5, where 
S v oc v a ), as measured between the 4.85 GHz Parkes-MIT-NRAO catalog (Griffith & Wright 1993) 
and the 1.4 GHz NRAO VLA Sky Survey (Condon et al. 1998, NVSS). This declination range is 
relatively unexplored for lenses and therefore likely to contain bright and useful specimens, and is 
far enough north to be within reach of the VLA and VLBA, instruments that facilitate the search. 

We obtained a snapshot of each object at 8.46 GHz with the VLA in its A configuration. 
Objects exhibiting multiple compact components (about 5% of the sample) were selected as lens 
candidates and scheduled for appropriate follow-up radio and optical observations. The goal of the 
follow-up observations is to determine whether the components have comparable spectral properties 
and surface brightnesses (as lensed images should have) and to search for a lens galaxy. 

This paper describes the second confirmed lens from our survey. The first was PMN J 1838- 
3427, described by Winn et al. (2000a). Following the organization of that paper, we end this 
introduction with a short chronology of the discovery of PMN J2004-1349, then examine the 
evidence in detail in subsequent sections. 

The initial VLA image of 1998 May 18 contained two compact components separated by 1'.'13. 
In the majority of such cases, higher-resolution radio imaging reveals one component to be flat- 
spectrum and compact, and the other to be steep-spectrum and extended: a core-jet structure. 
However, both components were unresolved in a MERLIN image obtained on 2000 April 1, which 
had 2.5 times the resolution of the VLA image. The components were also nearly unresolved in a 
VLBA image obtained on 2000 April 28, which had 100 times the resolution of the VLA image. 

With two compact components in such close proximity, PMN J2004-1349 was established as 
either a binary quasar or gravitational lens. 7 The most conclusive way to distinguish between these 
possibilities is to search for optical evidence of a lens galaxy. Such evidence was found in optical 
images obtained on 2000 July 25 with the Blanco 4m telescope at CTIO. 

Section 2 presents the radio properties of the system, and § 3 presents the optical properties. 



7 Here and elsewhere in this paper, we refer to the radio components as "quasars" because they are compact radio 
sources with pointlike optical counterparts, although strictly speaking they are not spectroscopically verified quasars. 



Within the latter section, § 3.1 presents the evidence for the lens galaxy, and § 3.2 provides photom- 
etry of PMN J2004-1349 and 5 nearby reference stars. In the last section we review the evidence 
that PMN J2004-1349 is a gravitational lens and use the observed properties to speculate on the 
nature of the lens galaxy and to estimate the order of magnitude of the time delay. 



2. Radio images 

We present results from radio images obtained with the VLA, MERLIN and VLBA, at fre- 
quencies ranging from 5 GHz to 22.4 GHz. In all the images, there are two components with nearly 
equal fluxes, separated by l'/13 at P.A. 60?6. We refer to the northeast component as NE and 
the southwest component as SW. Based on the average of our VLA measurements, the position of 
component NE is R.A. (J2000) = 20 h 04 m 07 s .091, Dec. (J2000) = -13°49'3lV07, within Qf.'l. 



2.1. VLA images 

We observed PMN J2004-1349 on 1998 May 18 for 35 seconds with the VLA at 8.46 GHz as 
part of the first phase of our lens search. Radio source 3C286 was used to set the absolute flux 
density scale, following the procedure suggested in the VLA Calibration Manual and adopting a 
flux density of 5.52 Jy. 

On 2000 November 1 we obtained additional VLA observations at 14.94 GHz and 22.46 GHz. 
These observations employed fast-switching between the target and the phase calibration source 
J2011-157 (from the VLA Calibration Manual). The total dwell time on PMN J2004-1349 was 
5 minutes at 14.94 GHz and 8.3 minutes at 22.46 GHz. To set the absolute flux density scale we 
observed J2355+498, a calibration source that is monitored monthly by S. Myers and G. Taylor of 
the NRAO and has proven to have a constant flux (within 5%) at these frequencies. The assumed 
flux densities of J2355+498 were 0.602 Jy at 14.94 GHz and 0.473 Jy at 22.46 GHz. 

In all cases, the total observing bandwidth was 100 MHz. Calibration was performed with 
AIPS using standard procedures, and imaging was performed with Difmap (Shepherd et al. 1994). 
In all cases, two compact components were obvious in the image prior to any deconvolution. We fit 
a model consisting of two elliptical Gaussian components to the visibility function using the "mod- 
elfit" utility of Difmap. After repeated iterations of model-fitting and phase-only self-calibration, 
we created a uniformly-weighted image with an elliptical Gaussian restoring beam. The 8.46 GHz 
image is shown in Figure 1. 

The component fluxes are listed in Table 1, along with the resolution and RMS level of the 
images. In every case, the separation between NE and SW was consistent with the more precise 
VLBA value listed in § 2.3. The Gaussian widths of NE and SW converged to very low values. 



2.2. MERLIN image 

On 2000 April 1 we observed PMN J2004-1349 with MERLIN, using all 6 telescopes of the 
array. The total observing bandwidth was 15 MHz centered at 4.994 GHz. To calibrate the antenna 
phases, we switched between the target and the nearby compact bright source J1939-1525 (from 
the VLBA Calibrator Survey) with a cycle time of 5 minutes. We divided the observation into 
four segments, separated in time by about an hour, in order to improve the ttu-coverage. The total 
dwell time on PMN J2004-1349 was 59 minutes. Radio source 3C286 was used to set the absolute 
flux density scale, assuming a flux density of 7.38 Jy on the shortest baseline. Calibration was 
performed at Jodrell Bank using standard MERLIN software, and imaging was performed with 
Difmap. 

Using the same procedure as was applied to the VLA data, we fit a two-component model 
to the visibility function. After repeated iterations of model fitting and phase-only self-calibration 
with a 5-minute solution interval, we arrived at our final image (Figure 1) and model (Table 1). 
The image is uniformly weighted and restored with an elliptical Gaussian beam. Both components 
were unresolved. 

The gravitational deflection of light is achromatic in the regime of lensing by galaxies, so the 
magnification ratio of point components should be independent of frequency. The near-equality 
of flux density ratios from 5 GHz to 22 GHz is therefore evidence that the two components are 
gravitationally lensed images. It is unlikely that two unrelated radio sources, or even the members 
of a binary quasar, would have identical continuum radio spectra. Stronger evidence for lensing is 
presented in § 3.1. 

We note that the flux density ratio measured at 8.5 GHz is not in formal agreement with the 
others, although it is within 2a. This discrepancy is not a serious challenge to the gravitational 
lensing interpretation. The 8.5 GHz measurement was made 2 years earlier than the others. Vari- 
ability of the background source causes the instantaneous flux density ratio to vary, because each 
lensed image displays the source at a different moment in its history. 



2.3. VLBA image 

We examined the system with all ten antennas of the VLBA on 2000 April 28. The phase 
center was approximately the midpoint between NE and SW. The observing bandwidth was 32 
MHz centered at 4.975 GHz. This bandwidth was divided into 4 intermediate frequencies, each of 
which was subdivided into sixteen 0.5 MHz channels. To calibrate the antenna phases, delays, and 
fringe rates, we switched between PMN J2004-1349 and J1939-1525 with a cycle time of 5 minutes. 
The total dwell time on PMN J2004-1349 was 66 minutes. 

Calibration, including fringe-fritting, was performed with standard AIPS procedures. We used 
a fringe-fit interval of 2 minutes on the phase reference source and applied interpolated solutions for 



delay and rate to PMN J2004-1349. Prior to imaging, we reduced the data volume by averaging in 
time into 6-second bins and in frequency into 1 MHz bins. This level of sampling was just enough 
to prevent significant bandwidth or time-average smearing over the required field of view. 

For imaging, we employed the multiple- field implementation of the Clean algorithm in AIPS, 
with a 512 x 512 field (with 20 mas pixels) centered on each component of PMN J2004-1349. 
The model constructed by the Clean algorithm was then used to self-calibrate the antenna phases 
with a solution interval of 2 minutes. This process, Cleaning and phase-only self-calibration, was 
iterated 10 times before arriving at the final model. The images based on this model, with uniform 
weighting and an elliptical Gaussian restoring beam, are displayed in Figure 1. The restoring beam 
is elongated in the north-south direction, with FWHM diameters 1.5 x 3.5 mas. 

Using the UVFIT task of AIPS, we fit a surface-brightness model consisting of two elliptical 
Gaussian components to the visibility data. To estimate of the uncertainty in each parameter, we 
observed the scatter in the results of separately fitting the model to individual 1 MHz spectral 
channels of data. The flux densities of the best-fit model components are listed in Table 1. Again, 
the flux density ratio is consistent with unity. The best-fit values of AR.A. and ADecl. are printed 
in the table caption. The models based on VLA and MERLIN data have separations that are 
consistent with these values, which are the most precise. 

The widths of component NE converged to values smaller than half the width of the restoring 
beam. We conclude NE is unresolved. Component SW converged to 2.06 x 0.76 mas at a position 
angle of —29°, with a scatter of 0.5 mas in each width and 16° in angle. Because the major axis is 
larger than the restoring beam, it appears that SW is resolved along one direction. These results 
are consistent with the visual appearance of the images in Figure 1. 

If component SW is truly resolved, then it has a lower surface brightness than NE, which 
complicates the hypothesis that PMN J2004-1349 is a gravitational lens. Gravitational deflection 
conserves surface brightness. Lensed images are therefore expected to have the same surface bright- 
ness in the absence of complications, such as propagation effects that are different for each image 
path from source to observer. 

The width of SW is not due to time-average or bandwidth smearing, which would in any case 
affect both components equally. Neither can the width be attributed solely to phase decor relation, 
because the width does not vary significantly as the self-calibration solution interval is altered, and 
because phase decorrelation would also affect both components. We conclude that SW is indeed 
resolved. This conclusion is worth checking with further VLBI observations, because the width is 
comparable to the resolution of the present data, and the uv -coverage of this short observation was 
poor. This point is discussed further in § 4.1. 



2.4. Total flux density 

Assessing the variability of gravitational lenses is important because flux monitoring of variable 
systems offers potentially significant scientific rewards, including the Hubble constant (Refsdal 
1964; Myers 1999; Schechter 2000) and the characterization of compact masses in the lens galaxy 
(Koopmans & de Bruyn 2000). Figure 2 is a logarithmic plot of the total flux density of PMN J2004- 
1349 as a function of radio frequency. It includes the entries from the PMN (Griffith et al. 1994) 
and NVSS (Condon et al. 1998) catalogs; the measurements with the VLA, VLBA and MERLIN 
described in this paper; and also two measurements with the Australia Telescope Compact Array 
(ATCA) that were taken on 2000 September 26 by J.E.J. Lovell (private communication). 

With the exception of the PMN flux density, the spectrum is a power law with a best-fit spectral 
index of a = —0.54 ± 0.10 (where S u ~ u a ). This is a typical value for radio-loud quasars. The 
large discrepancy between the PMN and the power law spectrum has two possible explanations. 
First, it could be the result of variability. The PMN flux is from 1990 whereas the others are 
from 2000. Second, the PMN flux is based on a single-dish measurement, whereas the others are 
interferometric. Emission that is smooth on scales larger than several arcseconds would be detected 
in the PMN survey but invisible to the interferometers. 

Variability is the more likely explanation. There is independent evidence that PMN J2004- 
1349 is variable on a time scale of years, from the discrepancy between the VLA (1998) and ATCA 
(2000) flux densities at 8.5 GHz. There is no independent evidence for extended emission nor is such 
emission expected from typical radio quasars. The NVSS survey used a 45" beam and might have 
been expected to detect this diffuse flux, but instead the NVSS flux lies along the same power-law 
line as the other interferometric measurements. There are no other sources in the NVSS catalog 
within 5' of PMN J2004-1349. 



3. Optical images 

On 2000 July 25 we obtained optical images of PMN J2004-1349 with the Blanco 4m telescope 
at CTIO. We used the Mosaic II CCD, which has eight 2048 x 4096 CCDs that are individually 
amplified and read out in pairs. We centered the target in chip #2. The night was photometric. 
Table § 2 reports the filter, exposure time, airmass, and seeing of each observation. 

The images from chip #2 were extracted and corrected for cross-talk from the paired amplifier. 
These images were then bias-subtracted and flat-fielded with standard IRAF 8 procedures. In addi- 
tion, for each of the two /-band images, sky fringes were removed by subtracting an appropriately- 
scaled fringe template for chip #2, which was kindly supplied by R.C. Dohm-Palmer. The two 



8 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



frames were then registered and added to create our final /-band image. The next section describes 
the optical counterpart of PMN J2004-1349 and, in particular, evidence of a lens galaxy in the 
/-band image. Section 3.2 provides photometry of PMN J2004-1349 and 5 nearby stars. 



3.1. Evidence of a lens galaxy 

The four panels of Figure 3 display the optical counterpart as viewed through each of the 
four filters. In I-band the quasars are easily distinguished, and it is apparent they have roughly 
the same separation, position angle and flux ratio as the radio double. In R the object is not as 
obviously double, and is not as symmetrical. The F-band counterpart is a single faint source lying 
between the expected quasar positions, which are marked with circles. In B the counterpart is 
barely detected if at all. 

Were PMN J2004-1349 a binary quasar, the optical counterpart would be double in all filters 
(with sufficient resolution and signal-to- noise ratio), and the separation would match the radio 
separation within the observational uncertainties. Our aim in this section is to demonstrate this is 
not the case. Rather, the actual morphology and its filter-dependence provides evidence of another 
light source between the quasars, which we believe is the lens galaxy. 

The most important evidence is in the /-band image, which is reproduced at higher contrast 
in Panel A of Figure 4. With the DAOPHOT package in IRAF, we constructed an empirical PSF 
of radius 6'.' 7, using the signal- weighted average of several bright and isolated field stars. We then 
fit a model consisting of two point sources to the optical counterpart of PMN J2004-1349. The 
best-fit separation was 880 ± 30 mas along a position angle of 62°, and the NE/SW flux ratio was 
1.01 ± 0.11. Compared to the radio double (1126.2 ± 0.2 mas, P.A. 61°) the optical double is along 
the same axis but is significantly smaller. 

Panel B is the residual image after the 2-point model has been subtracted. The positions of the 
subtracted components are marked with small circles. There are significant (4a) positive residuals 
just northwest and southeast of the axis of the double. A natural way to explain both the smaller 
optical separation and the pattern of residuals is that there is a diffuse source of unmodeled light 
centered between the quasars. The model components have been drawn together to account for 
this extra flux. 

To ensure that the diffuse flux is not an artifact of imperfect PSF fitting, we created an 
artificial double with the same position angle and separation as the radio double, using the image 
of an isolated star in the /-band image that is comparable in brightness to PMN J2004-1349. After 
applying the same PSF-fitting procedure as above, there was no trace of residual flux at the level 
seen in the actual data. 

Next, convinced of the reality of the diffuse flux between the quasars, we fit a two-point model 
to the original image in which the relative separation between the points was fixed at the VLBA 



value. The NE/SW flux ratio converged to 1.00 ± 0.11. Panel C of Figure 4 shows the residual 
image, in which the diffuse flux is seen to extend continuously across the axis joining the quasars. 
This diffuse flux is naturally interpreted as a lens galaxy. 

Based on the V and R morphology, the galaxy appears to be bluer than the quasars. This 
explains why there is a light source between the expected quasar positions in the F-band image, 
but the quasars themselves are not obvious. It may also explain why the i?-band components are 
not as clearly separate as in I, despite a comparable resolution and signal-to-noise ratio. This point 
is discussed further in § 4.2. 



3.2. Photometry of PMN J2004 1349 and field stars 

To represent the galaxy, we added a third point (Gl) to our model located exactly halfway 
between the quasars, and simultaneously solved for all three fluxes. There were still significant 
residuals to the north and south of Gl, so we added two additional points (G2 and G3) located at 
the peak residuals. Our final photometric model therefore contained two points representing the 
quasars (NE and SW), and three points for the galaxy (Gl-3). The relative separations were fixed 
at the values printed in Table 3, and the fluxes were allowed to vary simultaneously. Panel D of 
Figure 4 shows the I-band image with the quasars subtracted, in order to display the lens galaxy 
model. In Panel E, the galaxy model has been subtracted, highlighting the quasars. Finally, Panel 
F shows the residual image after all five components have been subtracted. We applied the same 
5-component model to the R- and y-band images. 

To establish a local photometric reference system for future observations, we used the same 
PSF template to find the instrumental magnitudes of 5 field stars. These stars are circled and 
labeled in Figure 5. For photometric calibration we observed the SA110 field described by Landolt 
(1992). Star #361 was used to determine the magnitude zero point, using an aperture diameter of 
14". We adopted "typical" CTIO extinction coefficients of ki = 0.06, kji = 0.11, ky = 0.15, and 
k B = 0.28 (Landolt 1992). 

The calibrated magnitudes of the reference stars and the components of PMN J2004-1349 are 
listed in Table 4. In V, the flux of component SE converged to zero, so we report only a lower limit 
on the magnitude. The non-detection in B implies B > 24.3 for all components. 

It is important to note that the quoted uncertainties are the statistical errors in the fitting 
procedure, and therefore apply only to magnitude differences. The calibration introduces an overall 
error of 0.05 magnitudes that affects all entries identically. Furthermore, the quoted uncertainties 
for the components of PMN J2004-1349 may be underestimated because they are internal to our 
choice of photometric model. Because all the components lie within two seeing discs, the magnitudes 
are covariant. 



4. Discussion 

4.1. Summary of evidence for lensing 

We now summarize the argument that PMN J2004-1349 is a gravitational lens. It consists 
of two radio components, each of which is compact on milliarcsecond scales and has a pointlike 
optical counterpart. This implies PMN J2004-1349 is a pair of radio-loud quasars, and invites 
three possible interpretations: a chance alignment of unrelated quasars, a binary quasar, or a pair 
of lensed images of a single quasar. 

The small separation (l'/13) argues against the chance alignment and binary quasar hypotheses. 
Using the analysis of faint 5 GHz radio source counts by Langston et al. (1990), the probability of 
such a close alignment of two unrelated sources brighter than 15 mJy is less than 2 x 10 -6 . The 
chance of finding such a pair among our initial sample of ~ 4000 radio sources is therefore less than 
1%. Furthermore, the small separation is more typical of lensed quasars (0 / /5 — 2") than binary 
quasars (> 3"; see e.g. Kochanek et al. (1999); Mortlock, Webster & Francis (1999)). 

The observation that the quasars have nearly the same flux ratio (~ 1:1) at 5 GHz, 8.5 GHz, 
15 GHz and 22 GHz would require a coincidence under the chance alignment and binary quasar 
hypotheses, but is a natural consequence of the lensing hypothesis. The strongest evidence for 
lensing is the presence of a diffuse source of light in the /-band image centered between the quasar 
positions, which is naturally interpreted as a lens galaxy. 

We believe this argument is conclusive. However, the lensing hypothesis is complicated by the 
observation that the quasars have different angular sizes in our VLBA images (§ 2.3). Because 
gravitational lensing conserves surface brightness, a pair of lensed images with equal flux should 
have equal angular sizes. Our data suggest that the SW quasar is resolved along one direction 
and the NE quasar is unresolved. If this observation withstands further scrutiny, there must be a 
propagation phenomenon that causes the surface brightnesses of the lensed images to differ. One 
possibility is differential scatter-broadening of the components by plasma in the lens galaxy or 
our own Galaxy. This hypothesis can be tested by determining whether the angular size varies 
characteristically as the square of the observing wavelength. One might also determine whether 
the percentage polarization or amount of Faraday rotation of the components are different, which 
would be independent evidence of different electron column densities. 



4.2. Nature of the lens galaxy 

The photometry described in § 3.2 was used to estimate the magnitudes of the lens galaxy 
(see Table 4). As described in that section, the quoted uncertainties are probably underestimates, 
because the galaxy is a low-surface-brightness object located between two quasars that are separated 
by only one seeing disc. It is even possible that what we have called the "lens galaxy" is more than 
one galaxy. Despite these caveats, our estimates can be used to speculate on the nature of the lens 
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galaxy, pending higher-resolution observations with adaptive optics or the Hubble Space Telescope. 
As a first step, we correct for Galactic extinction using the values of Schlegel et al. (1998), obtaining 
I = 21.5, R = 22.6 and V = 23.2 (within ~ 0.2 mag). 

Are these magnitudes consistent with the magnitude of a galaxy sufficiently massive to produce 
the observed image separation? The image separation produced by a lens galaxy with velocity 
dispersion a is roughly 

4 ^ = (1"4)( °- 

c 2 ; V220kms" 

The separation of l'/13 observed in PMN J2004-1349 therefore indicates a galaxy that is not much 

less massive than a typical L* galaxy. The magnitudes and colors of typical L* spiral and elliptical 

galaxies, over a range of redshifts, have been computed by Lehar et al. (2000) using several standard 

models for galaxy evolution. According to these calculations, an /-band magnitude of 21.5 implies 

a redshift between 0.5 and 1 (see Fig. 3 of Lehar et al. 2000). This is the range that is typical of lens 

galaxies, so the magnitude of the diffuse flux is certainly consistent with the lensing hypothesis. 



A0 ~ -^J- = (1?4) p r • (1) 



Furthermore, Lehar et al. (2000) show that for a given /-band magnitude, the V — I color is 
sometimes a useful diagnostic of galaxy morphology. In particular, for the redshift range 0.5 < z < 
1, early-type galaxies tend to be redder (V — I ~ 3) than late-type galaxies (V — I ~ 1.6). In this 
sense our lens galaxy photometry (V — I ~ 1.7 ± 0.3) favors a late-type galaxy. This statement 
depends upon the chosen model for the timing of star formation, and might be false for an early- 
type galaxies with a very late starburst. In addition, the conclusion may be sensitive to our choice 
of photometric model for the lens galaxy. Nevertheless, our working hypothesis is that the lens 
galaxy is a spiral galaxy with a redshift between 0.5 and 1. 



4.3. Magnification ratio and time delay 

In order to use a lensed quasar to determine the Hubble constant, one must not only measure 
the time delay At between the flux variations of the lensed images, but also construct a model 
of the gravitational potential of the lens galaxy that can be used to predict the quantity HoAt. 
The present observations do not allow for detailed mass modeling. Deeper and higher-resolution 
images, especially with adaptive optics or the Hubble Space Telescope, will be necessary for precise 
photometry and astrometry of the lens galaxy. Optical spectroscopy will also be essential, in order 
to obtain the redshifts of the background source and the lens. Despite the faintness of the quasars 
and the lens galaxy, low-resolution spectroscopy should not be a problem for an 8m-class telescope. 
In the meantime we can use the few observed properties of PMN J2004-1349 to predict the order 
of magnitude of the time delay. 

The simplest plausible model for a lens galaxy is a singular isothermal sphere, for which the 
time delay is 



A *=s(^ri< , + ; "< A » ,; 



fi + l 



(2) 
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where A9 is the image separation; [i is the magnification ratio; z\ is the lens redshift; and D[, D s 
and Di s are the angular-diameter distances to the lens, to the source, and between the lens and 
source, respectively. This expression can be easily derived from more general expressions (e.g., 
Witt, Mao & Keeton 2000). 

We further assume that the radio flux ratio, /, is a good approximation of the magnification 
ratio n, even though the flux ratio at any instant may depart from fj, due to variability in the 
background source, microlensing or interstellar scattering. A distinguishing feature of PMN J2004- 
1349 is that / ~ 1, which means we may only place an upper limit on the time delay. Judging 
from the values in Table 1, it seems likely that fi is within 5% of unity. For z s = 2 and z\ = 1 (the 
largest value consistent with the estimate in § 4.2), this implies hAt < 3 days. Smaller values of z\ 
produce smaller time delays. 

Such a short time delay would be challenging to measure, because it would require frequent 
sampling, and because it is poorly matched to the typical variability time scale of quasars (months). 
In any case, it may still be valuable to monitor this system to investigate microlensing within the 
lens galaxy, as has been done at both optical wavelengths (e.g. QSO 2237+0305; Irwin et al. 1989) 
and radio wavelengths (B1600+434; Koopmans & de Bruyn 2000). Here a small time delay is an 
advantage, because there is no need to wait a long time to obtain a light curve of the lagging image 
to compare with that of the leading image. 

As a final remark, we note that technically a singular isothermal sphere cannot produce a 
magnification ratio of unity. The symmetrical position, in which the lens is exactly aligned with 
the source, produces an Einstein ring. Small misalignments produce nearly equal doubles with high 
overall magnifications. If the magnification is indeed high, one might expect more sensitive radio 
images to reveal jets emerging from the quasar cores that form arcs or even an Einstein ring. These 
features would provide valuable constraints on lens models. 

Real galaxies have elongated potentials, which can also produce nearly equal doubles. The 
highly elongated potential of an edge-on spiral, in particular, tends to produce magnification ratios 
that are close to unity (Mailer, Flores & Primack 1997). This would be consistent with our suspicion 
based ovlV — I color (see § 4.2) that the lens galaxy is a spiral galaxy. In that case, lens modeling 
could provide interesting constraints on spiral galaxy structure, as it has for the edge-on spiral lens 
galaxy in B1600+434 (Koopmans, de Bruyn & Jackson 1998). 
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Fig. 1. — Radio images of PMN J2004-1349. In each panel, contours begin at 3a and increase by 
powers of 2, and the restoring beam is inset in the lower left of the image. The values of a and the 
beam dimensions are in Table 1. Upper left. VLA image (3"x 3"). Lower right. MERLIN image 
(3"x 3"). Lower left. VLBA image (20 mas x 20 mas) of component NE. Upper right. VLBA 
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Fig. 2. — Total flux density of PMN J2004-1349 as a function of radio frequency, on a log-log scale. 
The points labeled VLBA, MERLIN, ATCA and VLA are from measurements described in this 
paper. The points labeled PMN and NVSS are drawn from those radio catalogs. The solid line is 
the best fit after disregarding the PMN point, and has slope a = —0.54. Where error bars are not 
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Fig. 3. — Optical counterpart of PMN J2004-1349 as revealed by BVRI images from CTIO (see 
§ 3. Each panel is a 27" x 27" subraster with a logarithmic grayscale. North is up and east is left. 
In the V and B panels, the expected quasar locations are marked with circles. 
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Fig. 4. — Photometric model-fitting to the I-band image of PMN J2004-1349. Each panel is a 27" x 
27" subraster with a logarithmic grayscale. North is up and east is left. A. Original /-band image. 
B. Residual image after subtraction of a 2-component model in which the positions and magnitudes 
of the components were allowed to vary. The locations of the subtracted components are marked. 
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Fig. 5. — Wide-field (4' x 4') /-band image centered on PMN J2004-1349. The 5 reference stars 
discussed in § 3.2 are circled and numbered. Magnitudes and positions of these stars and of 
PMN J2004-1349 are printed in Table 4. The square box shows the field of view of the images of 
Figures 3 and 4. 



Table 1. Radio data for PMN J2004-1349 



Date 


Observatory 


Frequency 


Beam FWHM 


Flux 


density 


RMS noise 


Abs. flux 


Flux density 






(GHz) 


(mas x mas, PA.) 


N (mjy) 


S (mjy) 


(mjy /beam) 


uncertainty 


ratio 


1998 May 18 


VLA 


8.460 


307 x 167 (-20?7) 


14.41 


14.97 


0.26 


3% 


0.96 ±0.03 


2000 Apr 02 


MERLIN 


4.994 


147 x 37 (15?2) 


15.97 


14.52 


0.66 


5% 


1.10 ±0.07 


2000 Apr 28 


VLBA 


4.975 


3.5 x 1.5 (-0?2) 


16.37 


15.39 


0.22 


5% 


1.07 ±0.07 


2000 Nov 01 


VLA 


14.94 


176 x 108 (15?2) 


10.47 


10.06 


0.44 


5% 


1.04 ±0.06 


2000 Nov 01 


VLA 


24.46 


122 X 84 (22? 5) 


8.34 


8.07 


0.39 


10% 


1.03 ±0.09 



Note. — The most precise values for the separation between NE and SW come from the VLBA data: AR.A. 
: 552.34 ± 0.15 mas. All of the other radio data arc consistent with these values. 



981.47 ±0.08 mas, ADccl. 



CO 

I 



20 



Table 2. Journal of optical observations 



Filter 


Duration 


Seeing 


Airrnass 




(sec) 


(arcscc) 




1 


600 


1.07 


1.23 


I 


600 


1.15 


1.20 


B 


600 


1.13 


1.18 


V 


600 


1.00 


1.17 


R 


600 


1.24 


1.16 



Table 3. Photometric model 



Component 


AR.A. 


ADecl. 




(mas) 


(mas) 


NE 


490.7 


276.2 


SW 


-490.7 


-276.2 


Gl 








G2 


-28.5 


1028.7 


G3 


617.8 


-444.8 



Table 4. Photometry of PMN J2004-1349 and 5 reference stars 



Object 


AR.A. 


ADecl. 


I 


R 


V 




(arcscc) 


(arcscc) 


(mag) 


(mag) 


(mag) 


NE 


0.0 


0.0 


22.328( 92) 


23.639( 91) 


24.631(230) 


SW 


-1.0 


-0.6 


22.339( 84) 


23.929(155) 


> 24.6 


Gal. 


-0.5 


-0.3 


21.891(125) 


23.113(175) 


23.830(197) 


1 


108.6 


21.9 


19.289( 6) 


19.805( 6) 


20.432( 8) 


2 


27.2 


-43.6 


18.951( 5) 


19.738( 5) 


20.706( 7) 


3 


-0.2 


7.5 


19.912( 12) 


20.838( 9) 


21.901( 18) 


1 


-16.5 


78.8 


19.385( 7) 


20.994( 9) 


22.472( 25) 


5 


-72.9 


-12.5 


18.179( 7) 


19.002( 4) 


19.985( 5) 



Note. - - The parentheses contain the statistical error in millimagnitudcs, 
which does not include the overall calibration error of 0.05 magnitudes. The 
galaxy magnitudes are based on the sum of the fluxes of Gl, G2, and G3 
in the model of Table 3. Based on radio data, the coordinates of NE arc 
R.A. (J2000) = 20 h 04 m 07 s .091, Dec. (J2000) = -13°49'31'/07, within O'.'l. 



